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ABSTRACT 

We report the results of our observations of the S255IR area with the SMA at 1.3 mm in the very 
extended configuration and at 0.8 mm in the compact configuration as well as with the IRAM-30m 
at 0.8 mm. The best achieved angular resolution is about 0.4 arcsec. The dust continuum emission 
and several tens of molecular spectral lines are observed. The majority of the lines is detected only 
towards the S255IR-SMA1 clump, which represents a rotating structure (probably disk) around the 
young massive star. The achieved angular resolution is still insufficient for conclusions about Keplerian 
or non-Keplerian character of the rotation. The temperature of the molecular gas reaches 130-180 K. 
The size of the clump is about 500 AU. The clump is strongly fragmented as follows from the low 
beam filling factor. The mass of the hot gas is significantly lower than the mass of the central star. A 
strong DCN emission near the center of the hot core most probably indicates a presence of a relatively 
cold (< 80 K) and rather massive clump there. High velocity emission is observed in the CO line as 
well as in lines of high density tracers HCN, HCO+, CS and other molecules. The outflow morphology 
obtained from combination of the SMA and IRAM-30m data is significantly different from that derived 
from the SMA data alone. The CO emission detected with the SMA traces only one boundary of the 
outflow. The outflow is most probably driven by jet bow shocks created by episodic ejections from the 
center. We detected a dense high velocity clump associated apparently with one of the bow shocks. 
The outflow strongly affects the chemical composition of the surrounding medium. 

Subject headings: astrochemistry - HII regions - instrumentation: interferometers - ISM: clouds - 
ISM: molecules - stars: formation 


1. INTRODUCTION 

The formation of high mass stars (more massive than 
8-10 Mq) is still poorly understood. Studies of this pro¬ 
cess are currently a “hot topic” of astrophysical research. 
A major unsolved problem of high mass star formation 
is the characterization of accretion disks around young 
high-mass protostars. To date, most of evidenc e for the 
existence o f such disks has been indirect (e.g. Beuther 
et al. 2009). If present, these disks undoubtedly play a 
crucial role in the star formation process. Nevertheless, 
the very existence of these disks, much less their physi¬ 
cal properities, is not well-established. The importance 
of understanding these disks cannot be understated. For 
example, because of the larger masses involved, it is quite 
possible that high-mass disks are self-gravitating, un¬ 
like their Keplerian low-mass counterparts, implying that 
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there may be significant dynamical differences between 
the high-mass and low-mass star formation processes. 

S255IR is a part of the well known massive star form¬ 
ing complex located between the evolved Sharpless Hii 
regions S255 and S257. Several authors estimated the 
photometric distanc e to the complex at about 2.5 kpc 


(iRusseil et al. 2007 Chavarria et al. 2008 Ojha et al 
MTj r However, |Rygl et al.| ( |20I0[ ) report a distance of 


1.6 kpc based on trigonometric parallax measurements 
of methanol masers. This seems to be the most accurate 
distance estimate and we adopt it here, too. The previ¬ 
ous studies have sh own that S255IR contains a cluster o f 
early-B-type stars ( Howard et al. TT^IItoh e t al j200lD, 
several compact H ii regions dSnel l ^ bally |I98 d|), and 
a number of H 2 emission features (jlVIiralles et al. [1997 ). 
The total mass of the S255IR core has been estimated 
from single-dish observations at M ^ 300 — 400 Mq 
( Z inchenko et al.|[2009 Wang et ar]|2011 ). The estimate 
of the luminosity has been L ^ 2 x loT.© (|Wang et al. 


2011 ). 


The interferometric observations by Wang et al. (2011) 
revealed 3 continuum clumps and a high-velocity col¬ 
limated outflow in th e S255IR area. Our Paper I 
(Zinchenko et al.|2012) was devoted to the general struc- 
ture and kinematics of the complex based on the SMA 
observations at 1.3 and 1.1 mm ( the synthesized beam 
sizes were approximately 378 x 370 at 1.3 mm and 
279 X 276 at 1.1 mm) as well as VLA ammonia obser¬ 
vations and GMRT low f requency con t inuum data. Our 
results as well as data by Wang et al. (2011) indicate in 
particular a presence of a hot rotating core and a spectac¬ 
ular outflow in the S255IR region. However the angular 
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resolution was insufficient for a more detailed investiga¬ 
tion of this core. In addition, a reliable evaluation of 
physical and chemical properties of this system required 
observations of additional molecular transitions. There¬ 
fore we performed observations of the S255IR area with 
the SMA at a much higher angular resolution and at 
higher frequencies. 

Here we present observational data for S255IR ob¬ 
tained at 1.3 mm with the SMA in the very extended con¬ 
figuration and at 0.8 mm in the compact configuration. 
In addition we observed this area with the IRAM 30m ra¬ 
dio telescope in order to obtain short-spacing data com¬ 
plementing the SMA results. We mainly discuss proper¬ 
ties of the dense cores and high velocity outflows observed 
in this area. 


2. OBSERVATIONS AND DATA REDUCTION 
2.1. SMA 


The S255 IR area was observed with the Sub- 
millimeter Array (SMA) in its compact configura¬ 
tion on 2010 December 14th at 350 GHz. A three 
field mosaic was obtained with the following phase 
centers: 06^12^53L800, 17°59'22".l, 06^12^54L8876, 
17°59'29".3353 and 06^12^52L7124, 17°59'14".8647. 
The primary HPBW of the SMA antennas is 36" at these 
frequencies. Typical system temperatures on source were 
between 200 K and 400 K. The resulting uv coverage 
ranges from 10 kA to 90 kA (8 antennas in the array). 
3C454.3 and Uranus were used as the bandpass cali¬ 
brators and 0532+075, 0750+125 and 0530+135 were 
used as the complex phase and amplitude gain calibra¬ 
tors. A total of 8 GHz (342—346 GHz in the LSB and 
354—358 GHz in the USB) was observed with the SMA 
bandwidth doubling correlator configuration. The spec¬ 
tral resolution was 0.8125 MHz. 

In addition on 2011 January 07th S255IR was observed 
in the very extended configuration at 225 GHz. A sin¬ 
gle field was observed with the same phase center as for 
the central field at 350 GHz. The primary HPBW of 
the SMA antennas is 55" at these frequencies. Typi¬ 
cal system temperatures on source were between 90 K 
and 200 K. The resulting uv coverage ranges from 50 kA 
to 350 kA (7 antennas in the array). 3C454.3 and 3C279 
were used as the bandpass calibrators and 0750+125 and 
0530+135 were used as the complex phase and amplitude 
gain calibrators. A total of 8 GHz (217.0—220.8 GHz 
in the LSB and 229.1—232.95 GHz in the USB) was ob¬ 
served with the SMA bandwidth doubling correlator con¬ 
figuration. The spectral resolution was 0.406 MHz and 
1.625 MHz ( given the limited correlator capability, we 
were not able to have the high resolution across the whole 
band; so, we assigned different resolutions for different 
tracers). These data were combined with our previous 
data obtained at the same frequencies in the compact 
configuration (Paper I) which enables an investigation of 
the source structure in a wider range of spatial scales. 

The gain calibrator flux scale, calibrated against 
Uranus at 350 GHz and Callisto at 225 GHz, was found 
to be consistent within 5% with the SMA calibrator 
database and estimated to be accurate within 20%. 

The data calibration was carried out with the IDL su¬ 


perset MIR dScoville et al.||1993|), and subs equent imag- 
ing and analysis were done m MIRIAD dSault et al. 


ilysis 


1995). With robust weighting for the continuum and 
line data, the synthesized beam sizes are approximately 
272 X 179 at 350 GHz and 075 x 074 at 225 GHz. The rms 
noise is approximately 7 mJybeam”^ and 1 mJybeam”^ 
in the continuum images at 350 GHz and 225 GHz, re¬ 
spectively, and 100 mJybeam”^ and 20 mJybeam”^ in 
the spectral cubes at these frequencies at 2 kms“^ reso¬ 
lution. 


2.2. IRAM 30m radio telescope 

Single-dish observations of several molecular lines at 
the 30m IRAM radio telescope were performed in Oc¬ 
tober 2012 (N 2 H+ J = 3 — 2 at 279.5 GHz), December 
2012 (SiO J = 5 — 4 at 217.1 GHz) and January 2014 
(CO J = 3 — 2 at 345.8 GHz and CS J = 7 — 6 at 
342.9 GHz). The antenna HPBW was 973, 1179 and 775 
at these frequencies, respectively. The observations were 
performed in the OTF mode with position switching us¬ 
ing the HERA receiver at 217.1 GHz and EMIR receiver 
at higher frequencies. The central position of the maps 
was the same as for the primary field in the SMA ob¬ 
servations. The reference position was selected at —500" 
in right ascension from the central position. Apparently 
a weak CO(3-2) emission is present at the reference po¬ 
sition resulting in a weak negative feature in the CO 
spectra at Ulsr ~ 24 kms“^ (the bulk of the line emis¬ 
sion in this area is observed at Ulsr ^4 — 10 kms“^). 
This feature does not affect significantly the observed 
spectra. The system temperature was ^ 200 K for the 
N 2 H+ observations, ^ 300 K for the SiO observations 
and ^ 600 K for the CO and CS observations. Point¬ 
ing was checked regularly on nearby strong sources and 
pointing errors were within a few arcseconds. The an¬ 
tenna temperature calibration was made by the standard 
chopper-wheel method. 

The map size was approximately 2ffi x 2.'5 for the N 2 H+ 
and SiO observations, covering both S255IR and S255N 
regions, and approximately Iffi x l.'O for the CO and CS 
observations, covering only the S255IR area. Here we 
discuss only the data relevant to S255IR. The data on 
S255N are postponed for further publications. 

The data reduction was performed with the GILDAS 
package (http://www.iram.fr/IRAMER/GILDAS). 
Then the single-dish data were combined with the SMA 
data using the MIR IAD procedures as described by e.g. 
Wang et al. (2011). The conversion to the flux density 
scale was niade using the conversion factors presented 
on the IRAM 30m telescope website. 


3. OBSERVATIONAL RESULTS AND DATA ANALYSIS 

With the SMA we detected several tens of spectral lines 
in both 350 GHz and 225 GHz bands. A list of these 
lines, including their frequencies and energy of lower lev¬ 
els, is given in Tables The spectral line parameters 
are taken from the JPL ^ ickett et "ar]|1998[ ) and CDMS 
(jMiiller et al.||2001| |2005l) catalogs. 

We present the results in the form of maps as well as 
spectra and line parameters at selected positions. Eor 
continuum observations we give positions, flux densities, 
and size estimates of the continuum sources. 


3.1. Data analysis 
3.1.1. Methanol 
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Table 1 

List of molecular transitions observed at the SMA in 
the very extended configuration in S255IR. 


Table 2 

List of molecular transitions observed at the SMA in the 
compact configuration in S255IR. 


Molecule 

Transition 

Frequency 

El 



(GHz) 

(K) 


12C0 2-1 


2-1 

CH3OH 61 - 72 A- 
156 - I65 A- 
156 - I65 A+ 
2O1 — 2O0 E 
42 - 3i E 
80 - 7i E 
154 - I63 E 
8-1 — 7o E 
195 - 2O4 A+ 
195 — 2O4 A“ 
3-2 -4_i E 
222 - 21-3 E 
IO2 - 93 A- 
IO2 - 93 A+ 
I83 - 174 A+ 
CS 5-4 

DCN 3-2 

HNCO 10o,io — 9 o,9 

101.9 — 9i,8 

101.10 — 9i,9 
102,9 — 92,8 

102.8 — 92,7 

103.8 - 93,7 
103,7 — 93,6 

HC 3 N 24-23 

SO 65 - 54 

H 2 CO 3o,3 — 2o,2 

32,2 — 22,1 
32,1 — 22,0 

OCS 18-17 

19-18 

HCOOH 10o,io - 9 o ,9 

101.9 — 9i,8 

CH 3 CN 126 - lie 

125 - II5 
124 - II4 
123 - II3 
122 - II2 

12i - 111 

12 o - llo 


230.538000 

5.532 

220.398684 

5.289 

217.299202 

363.496 

217.642677 

735.160 

217.642678 

735.160 

217.886504 

497.919 

218.440063 

34.976 

220.078561 

86.051 

229.589056 

363.420 

229.758756 

78.076 

229.864121 

567.565 

229.939095 

567.561 

230.027047 

28.788 

230.292196 

598.499 

231.281110 

154.248 

232.418521 

154.248 

232.783446 

435.360 

231.220686 

22.194 

217.238530 

10.425 

219.798282 

47.471 

220.584762 

90.916 

218.981170 

90.569 

219.733850 

217.739 

219.737193 

217.739 

219.656770 

422.417 

219.656771 

422.417 

218.324711 

120.504 

219.949433 

24.429 

218.222195 

10.483 

218.475642 

57.608 

218.760071 

57.613 

218.903357 

89.304 

231.060983 

99.810 

220.038072 

48.061 

231.505705 

53.355 

220.594423 

315.313 

220.641084 

236.810 

220.679287 

172.555 

220.709017 

122.565 

220.730261 

86.849 

220.743011 

65.416 

220.747261 

58.272 


For the methanol data analysis, we constructed simple 
radiative transfer model, which uses the large velocity 
gradient (LVG) approximation. Dust emission and ab¬ 
sorption within the emission region was taken int o ac¬ 
count in the way described by Sutton et ah (2004). We 
assumed that the dust particles are intermixed with gas. 
The same physical temperature for the gas and dust com¬ 
ponents is assumed. The molecular emission region was 
assumed to be spherically symmetric and uniform in H 2 
density, gas and dust temperature, gas-to-dust ratio and 
methanol fractional abundance. The influence of exter¬ 
nal infrared sources was not considered. The dust opacity 
law was chosen as Tdust oc We adopted a gas-to- 

dust mass ratio of 100 and a cross sect ion at 1 mm of 
2.6 X 10“^^ cm^ (Sher wood et al.||l9 80). In addition to 
the model described in Sutton et al. "(21)04) we used the 
collision transition rates based on the model of collisions 
of methanol mole cules with He and para-H 2 molecules 
(Gragg et al.|2005 ). Scheme of energy levels in this model 
includes rotational levels with the quantum numbers J 
up to 22, \K\ up to 9; the levels include the rotational 
levels of the ground, first and second torsionally excited 


Molecule 

Transition 

Frequency 

(GHz) 

El 

(K) 


3-2 

345.795990 

16.596 

CS 

7-6 

342.882850 

49.372 

HCO+ 

4-3 

356.734242 

25.682 

CH 3 OH 

13i - 13o A-+ 

342.729796 

211.024 


13-1 - 14_2 E 

343.599019 

607.550 


I 82 - 173 E 

344.109039 

402.884 


19i - I 82 A+ 

344.443433 

434.697 


I 61 - 152 A- 

345.903916 

316.049 


18-3 - 17-4 E 

345.919260 

442.829 


13o - 12i A+ 

355.602945 

193.959 

HCN 

4-3 

354.505473 

25.251 


4-3 (0,li^,0) 

354.460435 

1049.892 


4-3 (0,li^,0) 

356.255568 

1050.021 


4-3 

345.339769 

24.861 

HC^^N 

4-3 

344.200109 

24.779 

HC 3 N 

38-37 

345.609016 

306.905 


39-38 

354.697456 

323.492 

SO 

83 — 77 

344.310792 

70.957 

33S0 

89,8 — 78,7 

343.086102 

61.564 


89,9 - 78,8 

343.087298 

61.566 

SO 2 

132,12 — 12 l,ll 

345.338538 

76.410 


124,8 — 123,9 

355.045517 

93.960 


134,10 — 133,11 

357.165390 

105.823 


154,12 — 153,13 

357.241193 

132.537 


114,8 — 113,9 

357.387580 

82.800 


84,4 — 83,5 

357.581449 

55.202 


94,6 — 93,7 

357.671821 

15.992 


74,4 — 73,5 

357.892442 

47.835 


64,2 — 63,3 

357.925848 

41.402 


174,14 — 173,15 

357.962905 

162.932 

H 2 CS 

10o,10 — 9 o,9 

342.946424 

74.132 


102,9 — 92,8 

343.322082 

126.830 


103,8 — 93,7 

343.409963 

192.613 


103,7 - 93,6 

343.414146 

192.556 

HCOOH 

16i,i6 — 15i,i5 

342.521225 

127.153 


15i,i4 — 14i,i3 

343.952413 

119.774 


16o,16 — 15o,15 

345.030596 

126.494 


162,15 — 152,14 

356.137250 

141.578 


states. In total, 861 levels of A-methanol a nd 852 levels 
of E-m ethanol were considered according to|Gragg et al.| 

([ 200 ^. 

With this model we made estimates of the hydro¬ 
gen number density, riHs, specific column density of 
methanol, A^cHsOh/AH, gas kinetic temperature, Tkin 
and fractional abundance of methanol A^ch 3 0 h/Ah 2 , us¬ 
ing the measured values of “quasi-thermal” (i.e. non¬ 
maser) methanol lines. The variable parameter of the 
source size is introduced in order to take into account 
beam-dilution effects. 

The brightness temperatures of all detected “quasi- 
thermal” methanol lines as well as upper limits for the 
brightness temperatures of other lines were taken into 
account. 

We have searched for a set of parameters which pro¬ 
vides the best agreement between the values of the calcu¬ 
lated brightness temperatures (TT^od^ and the measured 
brightness temperatures This corresponds to the 

minimum of ~ where is 

observational uncertainty for a particular line. 

The sources appear inhomogeneous in some cases. 
However, the LVG approximation treats the problem lo¬ 
cally and does not actually require physical homogeneity 
of the source. It can be well used for inhomogeneous 
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medium (e.g. Ossenkopf||1997). We have chosen the sim¬ 
plest approximation tor the case of the source which is 
smaller than the beam - introduced filling factor which is 
equal to the portion of the source size in the beam. The 
value of the filling factor is well constrained because the 
model is sensitive to its changes. We do not have enough 
observational constraints to study real dumpiness of the 
source. Introducing additional poorly constrained pa¬ 
rameters will reduce validity of the modelling. 


3.1.2. Methyl cyanide 

Methyl cyanide (CH 3 CN) is a symmetric-top molecule 
which is an efficient “thermo meter” for dense molecular 
gas (e.g. Boucher et al.|1980). To derive the kinetic tem¬ 
perature weTIseTthe'^opuIat ion diagram analysis which 
takes into account the optical depth in the CH 3 CN lines 
and th e beam filling factor as described by |Wang et al. 

(IMol). 


3.2. Millimeter wave eontinuum 


Maps of the continuum emission at 0.8 and 1.3 mm 
are presented in Fig. The 1.3 mm map was obtained 
by combining our new measurements in the very ex¬ 
tended configuration with the previous measurements in 
the compact configuration (Paper I). It shows a rather 
extended emission and compact cores in the SMAl and 
SMA 2 clumps. 

Estimates o f the continuum source parameters are 
given in Tables [3|^ For the measurements at 1.3 mm we 
indicate the parameters of the compact cores seen in the 
very extended configuration. Parameters of the more ex¬ 
tended components were presented in Paper I. The mea¬ 
sured fluxes are several hundreds mJy at 350 GHz and 
several tens mJy at 225 GHz in the very extended config¬ 
uration. The deconvolved sizes are from ^ 1 arcsecond 
to a few arcseconds at 345 GHz and about 0.3 arcsecond 
for the compact cores at 225 GHz. 


3.3. Basie properties of the moleeular emission 

Here we give the general description of the observed 
molecular emission. A more detailed information on the 
relevant species and transitions is presented in the fol¬ 
lowing sections. 

The general morphology and kinematics of the molec¬ 
ular emission was described in Paper I. Our new data set 
contains several tens of molecular transitions. In Fig. 
we present representative maps obtained with the SMA 
at 0.8 mm. 

These maps confirm that SMAl is the brightest source 
of molecular emission in this area. High-excitation lines 
are observed exclusively towards this clump. Several 
lower excitation molecular lines are detected also in 
SMA 2 . Emission of HCN, SO and CH 3 OH is observed 
also in SMA4. HGN and GS emission is probably present 
in the area to the west of SMA 2 , designated as S255IR- 
N 2 H+( 1 ) in Paper I. Since it is observed not in N 2 H+ 
only and for simplicity we shall designate it hereafter as 
SMA 2 -W. 

A new feature, not noticed in our Paper I and in other 
previous works is a molecular clump to the east of SMAl 
with a rather stroM emission in the HGN, HGO+, GS 
and SO lines (Fig. [^. We shall designate it as SMAl-E. 
This clump is located near the head of the jet observed 
in the NIR emission (see below). 


Spectra of several representative transitions tow ards 
the SMAl and SMA 2 clumps are presented in Figs. |3|4[ 
As mentioned in Paper I the main emission peak to¬ 
ward the SMAl clump is at Wsr ^4 — 5 kms“^. The 
line width is > 5 kms“^. Several lines including GO, 
HGN, HGO+, GS, SO show broad wings indicative of 
high velocity outflow. We discuss this feature in the fol¬ 
lowing sections. A rather strong emission is detected 
from vibrationally excited HGN with the excitation en¬ 
ergy ^ 1000 K above the ground level. 

The line emission from the SMAl-E clump peaks at 
about 8 kms“^. The linewidth is large, ^ 5 kms“^. 

The emission from the SMA 2 clump is observed at 
about 10 kms“^. The lines are narrow, ^ 2 kms“^. 

With the very extended array we detected the line 
emission almost exclusively from the SMAl compact 
core. No emission was detected in G^^O and SiO in¬ 
dicating an absence of compact structures in the lines 
of these molecules. The results of these observations are 
presented and discussed below. 


4. STRUCTURE, KINEMATICS AND PHYSICAL 
PROPERTIES OF DENSE CORES 

Our data presented in Paper I and here as well as 
data by Wang et al. (2011) indicate the presence of four 
continuum clumps in the S255IR area, designated from 
S255IR-SMA1 to S255IR-SMA4. In addition, in Paper I 
we detected clumps with a rather strong molecular emis¬ 
sion without continuum counterpart in the SMA data. 
One of them, S255IR-N2H+(1), is located close to SMA 2 
and was observed in the N 2 H+, NH 3 and several GH 3 OH 
lines. Estimates of their basic physical properties were 
presented in Paper I. Here we investigate further these 
objects using the new data set. 


4.1. SMAl 

SMAl is the brightest object in this area. As shown in 
Sect. [ 3 ] most of the observed molecular lines are detected 
only here. The deconvolved size of the continuum source 
at 350 GHz measured with the SMA in the compact con¬ 
figuration (Tableis close to that found in Paper I. The 
flux density measured at 350 GHz is only slightly higher 
than the flux density obtained at 284 GHz (Paper I). It 
is worth noting that in Paper I we could not separate the 
SMAl and SMA3 clumps. However, even if we take the 
integrated flux of these clumps at 350 GHz, the spectral 
index in the range 284 — 350 GHz will be only 1.6 which is 
inconsistent with a presumably optically thin dust emis¬ 
sion at these frequencies. Most probably this implies 
a significantly larger flux loss at 350 GHz in compari¬ 
son with the measurements at 284 GHz due to a smaller 
beam size. 

The deconvolved size of the continuum source detected 
with the SMA in the very extended configuration is about 
0'.'3 (Table which corresponds to about 500 AU. The 
observed aspect ratio for this core is close to unity. 


4.1.1. Kinematics 


Wang et al.| 

|( 

2011 ) noticed rotation of the core around 

the axis roughJ 

y parallel to the outflow direction. Now, 


at sub-arcsecond resolution we can investigate the core 
kinematics on smaller scales. In Fig. we present maps 
of the first moment of the CH 3 OH emission in the 42 — 3i 
E line and CH 3 CN emission in the I 23 — II 3 line. In this 
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Figure 1. Lef panel: map of the 0.8 mm continuum emission in the S255IR. The contour levels are (1, 2, 3, 4, 5, 6, 7, 8, 9, 
10)x40 mjybeam”^. The dashed rectangle indicates the area shown in the right panel. The dotted circles indicate positions of two 
molecular clumps discussed in the text. Right panel: map of the 1.3 mm continuum emission in the S255IR. The contour levels are (2, 3, 
4, 5, 7, 9, 11, 13) x4 mJv beam~^. The SMA beam is shown in the lower left corner of both panels. The crosses mark the position of the 
SMA3 clump according to |Wang et al.| ( [2011| ). 


Table 3 

Names, positions, flux densities, deconvolved angular sizes and 
position angles of the millimeter wave continuum sources measured 
at 350 GHz in the compact conflguration. 


Name 

a(2000) 

(h m s) 

5(2000) 

(0 , 

5350 

(Jy) 

^max 

(") 

^min 

(") 

P.A. 

n 

S255IR-SMA1 

6:12:54.00 

17:59:23.2 

0.50 

1.4 

0.6 

-11 

S255IR-SMA2 

6:12:53.76 

17:59:26.1 

0.56 

2.4 

2.1 

4 

S255IR-SMA3 

6:12:53.86 

17:59:23.7 

0.13 

1.3 

0.3 

44 

S255IR-SMA4 

6:12:54.01 

17:59:12.0 

0.25 

5.9 

3.8 

-82 


Table 4 

Names, positions, flux densities, deconvolved angular sizes and position 
angles of the millimeter wave continuum sources measured at 225 GHz 
in the very extended conflguration. 


Name 

a(2000) 

5(2000) 

‘S'225 

^max 

^min 

P.A. 


(h m s) 

(O , 

(Jy) 

(") 

(") 

(°) 

S2551R-SMA1 

6:12:54.010 

17:59:23.06 

0.058 

0.30 

0.27 

-19 

S2551R-SMA2 

6:12:53.779 

17:59:26.16 

0.016 

0.49 

0.31 

-27 


figure we also indicate the axis of the jet previously iden- 


tified through IR observations (Howard et ah 


l ocations of t he water masers measured with 
(Goddi et al. 


1997) and 


EeVLBA 


T'his hgure clearly shows that the core is really rotating 
around the axis of the jet. The rotation velocities along 
the line of sight amount to a few kms“^ (it varies from 
about 3 kms“^ to about 7 kms“^ for CH 3 OH 42 — 3i E 
and to about 8 kms“^ for CH 3 CN I 23 — II 3 ). In Fig.j^ 
we present the position-velocity diagram for the CH 3 OH 
emission in the 42 — 3i E line along the cut through the 
core center perpendicular to the jet axis (PA = 157°). It 
clearly shows the velocity gradient which can be consis¬ 


tent with Keplerian rotation but the angular resolution is 
still insufficient for firm conclusions about the character 
of this rotation. 

An implicit indication of a probable further increase of 
the rotation velocity in the innermost part of the core 
comes from the line width map (Eig. [^. It shows a sig¬ 
nificant increase of the line width in the center which 
can be related to rotation (the line width increases from 
about 2 kms“^ at the periphery to about 6 kms“^ in 
the center). 

It is also worth mentioning an increased line width 
along the jet axis (Eig. [^. It can be probably explained 
by an increased turbulence caused by the passage of the 
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06h2"'54f5 54f0 53f5 54f0 53f5 54f0 53f5 54f0 53f5 


a (2000) 

Figure 2. Maps of the integrated line emission for several representative transitions of various molecules (contours) overlaid on the image 
of the 0.8 mm continuum emission. The velocity interval is shown at the bottom of each panel. The dashed contours show negative features 
due to the missing flux. The contour levels are the following: (-3, -2, 2, 3, 5, 7, 9, 12, 15, 20) x A Jybeam“^ kms“^, where A = 3 for HCN 
and HCO+, A = 1 for H^^CN, HC^^N, HC 3 N, H 2 CS, CH 3 OH, SO 2 , and HCOOH, A = 2 for CS, and A = 1.5 for SO. The SMA 

beam is shown in the lower left corner of each panel. 


jet. 

The HCO+(4-3Wnd probably HCN(4-3) spectra to¬ 
wards SMAl (Fig.® show the red-shifted absorption dip 
which may be suggestive of infall. The CO(3-2) spec¬ 
trum measured with the SMA has a similar (although 
broader) feature (Fig.[^. At the same time the CO(3-2) 
spectrum obtained from the combined SMA and IRAM- 
30m data shows a peak at these velocities. Most proba¬ 
bly it means that this feature is related to an extended 
component resolved out by the SMA, although an infall 
cannot be fully excluded. 


4.1.2. Physical properties 

Let us consider physical parameters of this core. We 
estimated gas kinetic temperature from the CH 3 CN and 
CH 3 OH observations. Mod eling o f the CH 3 CN emission 
(Fig.[^ as described in Sect.|3.1 2 jyields kinetic tempera¬ 
ture in the range (74.9 - 197.6) K h cr confidence interval) 


towards the emission peak (TableWe have detected a 
large number (15) of CH 3 OH transitions towa rds SM Al 
(Table [T|. From the CH 3 OH analysis (Sect. 3.1.1) we 
obtain ^in = 178 K with the 165 — 195 K la confidence 
interval at the center, Tkin = 170 — 200 K at 072 to the 


Table 5 

Kinetic temperature, CH 3 CN column density, 
beam Ailing factor and line optical depth derived 
from the CH 3 CN data obtained in the very 
extended conflguration (HPBW O'/A) at the 
central position in SMAl. The Icr confldence 
intervals are indicated in the parentheses. 


Parameter 

Fit results 

Tkin (K) 

128.8 (74.9 - 

197.6) 

ArL(CH3CN) (IQi® cm- 2 ) 

2.12 (0.96 - 

5.80) 

Beam filling factor 

0.23 (0.18 - 

0.36) 

t{K = 0) 

2.53 (0.93 - 

8.37) 


north and Tkin = 140 — 165 K at 072 to the south (Ta¬ 
ble . Therefore there is probably a temperature gradi¬ 
ent m the north-south direction. The temperature seems 
to be anti-correlated with the CH 3 OH column density 
(Table®. 

We aaopt the gas kinetic temperature of 170 K as the 
average value (it is close to the weighted average of the 
CH 3 OH and CH 3 CN results). Assuming the same tem¬ 
perature (170 K) for the dust we obtain a total mass 
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Vlsr (km/s) 

Figure 3. Spectra of several representative molecular transitions towards the SMAl clump. The CO(3-2) spectrum is scaled by the factor 
of 0.5. It is obtained from the combined SMA and IRAM 30m data. The other spectra are from the SMA data only. The angular resolution 
is about 2 arcseconds. 
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Figure 4. Spectra of several representative molecular transitions towards the SMA2 clump. The CO(3-2) spectrum is scaled by the factor 
of 0.5. It is obtained from the combined SMA and IRAM 30m data. The other spectra are from the SMA data only. The angular resolution 
is about 2 arcseconds. 


Table 6 

Physical parameters derived from the CH 3 OH data obtained in the very extended configuration 
(HPBW OO 4 ) at several positions in SMAI. The Icr confidence intervals are indicat ed i n the 
parentheses (the density is not constrained by the model as evident from Fig. |I0|). 


Offsets 
(",") 

'Tkin 

(K) 

A^chqOh/AR 
(10^^ cm“^s) 

(10® cm-®) 

beam filling factor 
(%) 


0, +0.2 

183 (170—200) 

3.6 ( 2 . 8 — 5 . 0 ) 

3.2 

14.8 (13.9—16.4) 

10-8(> 10“'^) 

0, 0 

178 (165—195) 

5.6 (4-0—9.5) 

3.2 

16.0 (15.2—17.2) 

10-8(> lO-’^) 

0, -0.2 

153 (140—165) 

8.9 (5.0—12.6) 

0.2 

15.2 (13.3—16.0) 

10-8(> 10-7) 


of this hot component of about 0.3 Mq. As in Paper I 
we assume a gas-to-dust mass ratio of 100 , and adopt 
a dus t abso rption efficiency following [Ossenkopf fc Hem] 


mug 


(1994). 


The peak gas column density estimated 
from the continuum data is N{ll 2 ) ^ 3 x 10^^ cm“^. 
With a size of 500 AU the mean density of hot gas is 
about 6 X 10^ cm“^. This is 2 times higher than our es¬ 
timate of the SMAl mean density in Paper I taking into 
account the difference in the adopted distances here and 
in Paper I. This increase can be expected since now we 
consider much smaller scales near the core center. The 
CH 3 OH modeling puts no significant constraints on den¬ 


sity. The relative abundance of CH 3 OH is ^ 10“^. 

Both CH 3 CN and CH 3 OH observations indicate beam 
filling factor of 0.15-0.2. This means that the source 
is very inhomogeneous on the 074 scale and probably 
consists of clumps with the size <C 500 AU and density 
much higher than the mean density found above. This 
higher density estimate does not contradict the CH 3 OH 
data. A similar picture of clumpy medium was inferred 
from our obser vations of high mass star forming regions 
on lar ger scales (Pirogov & Zinchenko 2008 Pirogov et al. 
2012 |). Most probably it reflects turbulence in the cores 


Uur data set includes several other tracers of the hot 
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Figure 5. Maps of the first moment of the CH 3 OH emission in the 42 — 3i E line (left panel, color scale) and CH 3 CN emission in 
the 123 — II 3 line (right panel, color scale) towards the SMAl clump. The intensity cut-off is 100 mjybeam“^ kms~ ^ for CH 3 OH and 
200 mJybeam“^ kms“^ for CH 3 CN. Contours show the 1.3 mm continuum emission. The triangles mark water masers (iGoddi et al.|2007t 
and the arro ws indicate their proper motions (in the cases when they are measured). The red line indicates the jet axis as tound by|llowai^| 
|et al.| (|1997|). The beam for the molecular maps is shown in the lower left corner of both panels. The scale for the velocities of the proper 
motions IS shown in the lower right corner. The pixels with the lowest velocities appear as the blue squares in the right panel. 



Figure 6. The position-velocity diagram for the CH 3 OH emis¬ 
sion in the 42 — 3i E line along the cut through the SMAl center 
perpendicular to the jet axis (PA = 157°). The dashed curves 
correspond to Keplerian rotation around a central mass of 20 Mq 
with the inclination angle of 25°, the dotted curves correspond to 
the inclination angles of 25° ± 5°. 



a(2000) 

Figure 7. Map of the second moment of the CH 3 OH emission 
in the 42 — 3i E line (color scale). Contours show the 1.3 mm 
continuum emission. The beam for the molecular map is shown in 
the lower left corner. 


gas. One of the most important is HNCO. We detected 
HNCO lines in different K-i ladders (LC-i = 0,1,2^) 
with the excitation energies up to ^ 400 K (Table [B. 
In Fig. we present the rotational diagram for the ob¬ 
served HNCO transitions. This diagram was obtained 
from peak integrated line intensities found by 2D Gaus¬ 
sian fitting of the integrated intensity maps in different 
HNCO lines. We used the line strengths and dipole mo¬ 
ment fro m the Cologne Database for Molecular Spec¬ 
troscopy (jMuller et al.||2001l |2005|). 

This diagram indicates a rotational temperature of 
318 ± 70 K (Table [^. However, the K-i = 0,1 transi¬ 
tions can be saturated. The peak brightness temperature 
in these transitions is close to the brightness tempera¬ 


ture in the apparently optically thick lines (e.g. CO and 
CHsCN) an d our modeling using RADEX (|van der Tak 
et al.|[2QQ7 ) shows that at the derived physical parame- 


ters and column density the optical depth in these lines 
is about unity. In this case the derived rotational tem¬ 
perature would apparently represent an upper limit to 
the excitation temperature. The deconvolved source size 
is - 0'.'3 X 072 for the K-i = 0 transition and decreases 
for higher K-i ladders. The map of the integrated in¬ 
tensity in the HNCO K_i = 2 transition is shown in 
Fig, along with the map of the OCS emission which is 
another tracer of hot gas. The peaks of the HNCO and 
OCS emission practically coincide with the continuum 
peak, although distributions of these molecules seem to 
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-20 0 20 40 

Vlsr (km s ) 

Figure 8 . The SMA spectra of HCO+(4-3) and CO(3-2) as well 
as the combined SMA+30m spectrum of CO(3-2) towards SMAl. 
The angular resolution is about 2 arcseconds. 



Figure 9. The CH3CN spectrum measured towards the SMAl 
emission peak and the model fit. 


Table 7 

Fit results for the rotational diagrams 
obtained towards SMAl. The Icr 
uncertainties are indicated in the 
parentheses. In the last column the 
approximate HPBW for the corresponding 
data is indicated. 


Molecule 

Tkin 

(K) 

log Al 
( cm“^) 

HPBW 

(") 

HNCO 

318(70) 

16.34(0.04) 

^OA 

SO 2 

146(16) 

15.50(0.05) 

^2 


65(11) 

15.35(0.08) 

^2 


be somewhat different. 

The total HNCO column density derived from the ro¬ 
tational diagram is ^ 2 x 10^^ cm“^. We used the par¬ 
tition function from the Cologne Datab ase for Molecular 
Spectroscopy ( Muller et al.||2QQl 2005), too. A compar- 



160 170 180 190 200 

Tk K 


Frequency MHz 


Figure 10. An example of the CH 3 OH data fitting for the central 
position in SMAl. Variations of with the model parameters are 
shown in the 2 top and left bottom panels. The minimum of x^ 
is marked by the star. Confidence intervals of 0.25, 1, 1.6 2.6 
cr are plotted on the grey scale. The confidence interval of 1 cr is 
marked by the white line. The comparison of the model (stars) and 
observed main-beam temperatures is shown in the right bottom 
panel. 



100 


200 

Eu (K) 


300 


400 


Figure 11. The rotational diagram for HNCO towards the emis¬ 
sion peak. 


ison with the total gas column density indicates HNCO 
abundance of X(HNCO) ^ 10“^. This is a rather high 
value, close to the highest HN CO abundance derived in 
the survey of massive cores by Zinchenko et al. (2000). 


However, taking into account tne note above about the 
HNCO lines saturation, this value may need a correction. 

In the 350 GHz band we detected many SO 2 lines (Ta¬ 
ble!^. The corresponding population diagram is shown 
in Fig. It indicates a range of temperatures. For 
low-excitation transitions the rotational temperature is 
65 ± 11 K. For the high-excitation ones it is 146 ± 16 K. 
The SO 2 observations were performed at a much lower 
angular resolution (2 arcsec) than in case of HNCO and 
apparently include emission from both the hot core and 
the surrounding cloud. The derived temperatures are 
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Figure 12 . The maps of the integrated intensity of HNCO Jk_i = IO 2 — 82 (left panel) and OCS J = 18 — 17 (right panel) transitions 
(color scale) towards the SMAl clump. The intensity units are Jybeam“^ kms“^. The beam size is 074. The contours show the 1.3 mm 
continuum emission with the same angular resolution. The beam for the molecular maps is shown in the lower left corner of both panels. 

tically the same as those of high excitation lines of other 
molecules tracing apparently the hot gas. Therefore, 
most probably the observed DCN emission arises within 
the disk and since we probably see the disk nearly face- 
on as discussed above, the physical distance from the 
center is not much larger than the projected distance. 
An estimate of the abundances in the LTE approxima¬ 
tion gives X(DCN) > 10“^^ and X(^^CS) > 3 x 10“^^ 
(assuming the total mass of the clump < 10 Mq and 
temperature of 50-100 K). This is a “normal” value for 
while the derived DCN abundance implies a signifi¬ 
cant deuteriation as follows from compar ison with typical 
HCN abundances in massive cores (e.g. [Zinchenko et al 


' ' 

' 1 ^ 

^ 1 1 ^ ^ 1 1 1 ^ = 






. ' 


Eu (K) 

Figure 13. The population diagram for SO 2 at the SMAl emis¬ 
sion peak. 

consistent with the other estimates for these components. 

Another indicator of the hot environment is the vibra- 
tionally excited HCN. We detected the emission in the 
i ;2 = 1 state, about 1000 K above the ground state. For 
the V 2 = 2 emission the upper limit is about 5 times 
lower. Follo wing the analysis presented in e.g. IVeachj 


et al. (2013) we obtain an upper limit for the excitation 


temperature between these states of about 500 K. This is 
consistent with the other estimates of the gas tempera¬ 
ture given above. The critical density for HCN excitation 
is high (> 10^^ cm“^, Veach et al.|2013 ) which is consis¬ 
tent with our density estimates tor the hot gas. 


2009|). As poin t ed ou t in Paper I, recent modelling byjAl 


bertsson et al. (2013) shows that DCN/HCN abundance 
ratio sharply drops at temperatures > 80 K. It means 
that temperature of the DCN emitting clump should be 
rather low, much lower than the temperature of the hot 
gas in the disk. Another possible explanation, as also dis¬ 
cussed in Paper I, can be a very young age of the clump, 
insufficient to reach the steady-state DCN/HCN abun¬ 
dance ratio, but this looks less probable. In principle, 
interiors of the accretion disk can be well shielded from 
an external radiation, probably providing necessary low 
temperatures. In any case this clump deserves further in¬ 
vestigation in order to clarify its properties and nature. 
It is worth mentioning that the mass of this clump cannot 
be lower than ^ 1 Mq, otherwise the inferred ^^CS abun¬ 
dance would be unrealistically high. The virial mass we 
estimate about the same, hence the clump can be gravi¬ 
tationally bound. 


4 . 1 . 3 . A cold clump in the hot core? 


In Fig. M we plot maps of the DCN J = 3 —2 and ^^CS 
J = 5 —4 integrated line emission in the SMAl core. The 
emission regions are very compact and the emission peak 
is shifted from the continuum peak which apparently co¬ 
incides with the YSO location. The projected distance 
from the continuum peak is roughly 300 AU. 

The velocities of the DCN and ^^CS emission are prac¬ 


4 . 2 . SMA2 

As in case of the SMAl the spectral slope for the con¬ 
tinuum emission in the frequency range 284-350 GHz is 
too low 1.3), which is probably caused by a higher 
flux loss at higher frequencies. With the SMA at sub- 
arcsecond resolution we detected a compact component 
in continuum emission of about the same size as in 
the SMAl. However, no high-excitation molecular lines 
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Figure 14. The maps of the integrated intensity of the DCN J = 3 —2 (left panel) and J = 5 —4 (right panel) transitions (color scale) 

towards the SMAl clump. The intensity units are Jybeam“^ kms“^. The beam size is O'.'d. The contours show the 1.3 mm continuum 
emission with the same angular resolution. The beam for the molecular maps is shown in the lower left corner of both panels. 

could be detected in this area. In Paper I we derived the 
temperature of the SMA 2 clump ^ 40 K. Now with the 
exten ded da ta set, from the methanol excitation analysis 
(Sect. |3. 1 . 1 ) we found the la confidence interval for tem- 
perature of 40-80 K. The best estimate of the methanol 
relative abundance is ^(CHsOH) ^ 10“^. These esti¬ 
mates are obtained at the 2 arcsec scale. It cannot be 
excluded that the temperature of the compact structure 
observed in the very extended configuration is somewhat 
higher. However it cannot be much higher since there is 
no sign of a higher temperature component in the molec¬ 
ular data. Assuming the temperature of 50 K we obtain 
the mass of the compact core ^ 0.2 Mq. 

To the west from the SMA 2 there is an area of molecu¬ 
lar emission without continuum counterpart in the SMA 
data, observed earlier in the N 2 H^, several 

CH 3 OH lines (Paper I andjWang et al.||2011|). Our new 
methanol data analysis indicates the temperature in the 
range 25-65 K, density n(H 2 ) in the range 6.3 x 10^— 

2.5 X 10^, ^(CHsOH) ^ 10“^. Some of the methanol 
transitions are inverted in the model and can be masing. 

These high gas densities and rather high methanol abun¬ 
dances at relatively low temperatures can be explained 
by the influence of shock. 


4.3. SMA3 

In Paper I we did not detect this compone nt which 


coincides wi th the near-infrared source NIRS 1 (Tamura 
et al. 1991) and is identified as a massive disk candi- 


(Fig-i or molecular emission at this position. 

4.4. SMA4 

This clump was first detected in Paper I. It shows 
a weak continuum and spectral line emission in sev¬ 
eral molecular transitions. The measured continuum 
flux density at 350 GHz is practically the same as at 
284 GHz which also implies a much higher flux loss at 
the higher frequency. The temperature derived from our 
new methanol data analysis is about 45 K, higher than 
estimated in Paper I from ammonia data. However the 
ammonia emission here is very weak and the uncertain¬ 
ties are high. Some methanol transitions are inverted in 
the model. The influence of shock is also probable. 

5. MORPHOLOGY AND PROPERTIES OF THE OUTFLOWS 

In the S255IR area high velocity emission is detected 
in lines of GO, SiO and several high d ensity tracers in¬ 
cluding HGN, HGO+ and GS. In Figs. 15|l6 we present 
channel maps of the GO J = 3 — 2 emission in the blue 
and red line wings, respectively. The maps of the in¬ 
tegrated line wing emission (Fig. ^7) show that the GO 
high velocity emission observed witTithe SMA looks like a 
highly collimated bipolar outflow originating near SMA 2 . 
There is another more compact component near SMAl. 

On the other hand the bipolar outflow observed in 
HGN and HGO+ is apparently associated with SMAl 
while is parallel to the GO flow (Fig. pT|). It is appar- 


date by N IR polarization observ ations (Jiang et al. 2008) 


However, Simpson et al. (2009) argued that these mea¬ 
surements could be attected by instrumental effects and 
there is no real evidence for the “polarization disk” here. 
At the same time the NIR measurements indicate a pos¬ 
sible outflow related to this object. A c omparison with 
the measurements by Wang et al. (2011) gives a spectral 
index of about 3.3 in the frequency range 225-35 0 GHz, 
consistent w ith the optically thin dust emission. |Wang| 
et al. (2011) estimated mass of this clump at about 2 Mq. 
At sub-arcsecond resolution we do not see any continuum 


ently associated with the jet observecTm particular in 
the Fell emission which is also plotted in this figure. It 
is worth mentioning that the extent of the jet is much 
smaller than the extent of the outflow seen in CO, al¬ 
though their orientations coincide. 

The SMA data alone hint at two parallel outflows with 
different origins. However single-dish CO(3-2) observa¬ 
tions with the IRAM 30m telescope show a different pic¬ 
ture. Here we see a wider outflow clearly originating at 
SMAl (Fig. While the emission peaks coincide in 
the SMA ano^m maps, near the driving source the CO 
emission measured with the SMA traces only the north- 
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Table 8 

Physical parameters derived from the CH 3 OH data obtained in the compact configuration (HPBW 3"). The Icr 

confidence intervals are indicated in the parentheses. 


Object 

Tkin 

(K) 

-^cHoOh/^V' 

(10i^cm-3s) 

nu2 

(cm-3) 

beam filling factor 
(%) 

A^CHaOn/A^Hs 

SMA 2 

SMA2-W 

SMA4 

50 (40—80) 
40 (25—65) 
45 (30—90) 

2.0 (1.8—3.5) 
3.2 (2.2—4.5) 
3.6 (0.8—7.9) 

6.3x10® (< 10®) 

1.0x10® (6.3 X 10^—2.5 X 10®) 
5.6x10^ (103—6.3 X 10®) 

99.3 (> 84) 

64.0 (63—70) 

99.9 (> 20) 

10 -®(> 10“®) 
10 -'^(> 10“®) 
10 -'^(10-®—lO-'^) 



6^12'^56^ 55" 54" 53" 52" 51" 55" 54" 53" 52" 51" 55" 54" 53" 52" 51" 55" 54" 53" 52" 51" 


RA (J2000) RA (J2000) RA (J2000) RA (J2000) 

Figure 15. Channel maps of the CO blue wing emission in the S255IR. The numbers in the upper left corner indicate the channel 
velocity in kms“^. The dashed contours show negative features due to the missing fiux. The contour levels are (-3, 3, 5, 7, 10, 15, 20, 30, 
40) X100 mJybeam”^. The SMA beam is shown in the lower left corner of each panel. 
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Figure 16. Channel maps of the CO red wing emission in the S255IR. The numbers in the upper left corner indicate the channel velocity 
in kms“^. The dashed contours show negative features due to the missing fiux. The contour levels are (-3, 3, 5, 7, 10, 15, 20, 30, 
40) X100 mJybeam”^. The SMA beam is shown in the lower left corner of each panel. 
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Figure 17. Left panel: maps of the CO(3-2) high velocity emission as observed with the SMA (blue and red contours) in the S255IR 
area overlaid on the continuum image at 0.8 mm. The SMAl and SMA2 continuum clumps are marked. The velocity intervals are from 
-63 kms“^ to -13 kms“^ for the blue wing and from 25 kms“^ to 75 kms“^ for the red wing. The dashed rectangle indicates the area 
shown in the right panel. Right panel: maps of the HCO+(4-3) high velocity emission (blue and red contours) in the S255IR area overlaid 
on the continuum image at 0.8 mm. The velocity intervals are from -1 9 kms~^ to 3 km s~^ for the blue wing and from 11 kms“^ to 
27 kms“^ for the red wing. The black contours show the Fell emission |Wang et al.|201l|. The map of the 15 GHz continuum emission 
(from the VLA archival data, the angular resolution is about 4 arcsec) is plotted wifli green contours. The dashed line indicates the jet 
axis (PA = 67°) as in Fig.[^ The SMA beam is shown in the lower left corner of both panels. 
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a (2000) 

Figure 18. Maps of the CO(3-2) high velocity emission as ob¬ 
served with the IRAM 30m telescope (blue and red thick contours) 
in the S255IR area overlaid on the continuum i mage at 0.8 mm. 
The thin contours show the SMA maps from Fig. |17| The dashed 
line indicates the axis of the position-velocity cut (PA = 67°). The 
orange contours show the N2H+(3-2) integrated line emission ob¬ 
tained by combining the SMA and 30m data. T he c ross marks 
the position of the high velocity dense clump (Sect. [5l^ . The 30m 
beam and the SMA beams for the CO and N 2 H+ observations, 
respectively, are shown in the lower left corner (from left to right). 


ern edge of the outflow. This shows that the SMA map 
gives a distorted picture aud iu reality we have here a 
wider augle outflow origiuatiug at the SMAl. 

Appareutly this meaus that this uortheru edge of the 
outflow coutaius a relatively bright compoueut with a 
characteristic scale comparable to the SMA beam. The 
southeru edge of the outflow uear the driviug source 
would be more diffuse aud resolved out by the SMA. Most 
probably this implies that the CO emissiou is formed iu 
a compressed layer surrouudiug the outflow cavity. For 
some (uuclear) reasou iu the uortheru part this layer is 
more prououuced. It is worth meutiouiug that the bright¬ 
est peaks of the high velocity CO emissiou practically 
coiucide iu the SMA aud 30m maps aud also coiucide 


with spots of molecular hydrogeu emissiou (Waug et al 


2011) which probably iudicate bow shocks. Basic phyi 


ical pyame ters of the outflow were estimated by |Waug| 


et al. (2011) from their CO( 2 -l) data. Our uew data are 


cousisteut with these estimates. 

The positiou-velocity (P-V) diagram for the outflow 
aloug the cut iudicated iu Fig. coustructed from the 
30m data, is showu iu the left pauel of Fig. 19 The gap 
at IAsr ^ 24 kms“^ is caused by the emissiou at the 
refereuce offset positiou meutioued above. Iu the right 
pauel we preseut the P-V diagram obtaiued from the 
SMA CO(3-2) data. The maiu features of these diagrams 
coiucide. They show uou-mouotouic depeudeuce of the 
velocity ou offset, which hiuts at two outflow compoueuts 
at differeut distauces from the driviug source. These P-V 
diagrams are some what diffe reut from that prese uted by 


Waug et al. (2011). However Waug et al. (2011) plotted 
this diagram tor a differeut cut, with the positiou augle 
of 75°, which makes the direct comparisou irrelevaut. 


Iu Fig. 20 we plot the P-V diagram for the HCO+(4- 


3) emissiou as observed with the SMA, aloug with the 
part of the CO(3-2) P-V diagram for the same iutervals 


of the offset aud velocity. The diagrams are very similar, 
which meaus that both molecules trace appareutly the 
same gas. 


5.1. Dense high veloeity elump 

Our data show a stroug, compact blue-shifted CS aud 
HCN emissiou close to the peak of the CO blue-shifted 
hue wiug emissiou. The peak of this CS aud HCN emis¬ 
siou is marked by cross iu Fig.f^ The CO(3-2), CO( 2 - 

1 ) , HCN(4-3), CS(7-6), HC(T^(4-3) aud N 2 H+^ 2 ) 
spectra towards this positiou are preseuted iu Fig. 

Oue cau see rather stroug aud broad HCN aud CSTInes 
at a ceutral velocity of about -5 km s“^. There is a hiut of 
wiugs iu these hues exteudiug to about -20 kms“^. The 
peak of the CO emissiou, especially iu the higher J = 
3 — 2 trausitiou is observed at more uegative velocities. 
There is uo detectable HCO+(4-3), N 2 H+( 3 - 2 ), CH 3 CN 
aud coutiuuum emissiou (for N 2 H+ aud CH 3 CN we use 
the data from Paper I). 

The spatial distributiou of the CS(7-6) aud CO(3- 

2 ) emissiou iutegrated iu the velocity rauges - 
13...+1 kms“^ aud -19...+1 kms“^, respectively, is 
showu iu Fig. The HCN distributiou is very simi¬ 
lar to that of CST The decouvolved size of the CS emit- 
tiug clump is about I'.'l x 073 which correspouds to 
1800 AU X 500 AU. It is appareutly located at the head 
of the stream observed iu CO, almost exactly at the jet 
axis showu iu Fig. 

At the velocity ofthe CS aud HCN emissiou peak, the 
brightuess temperatures iu the CO(3-2), CO( 2 - 1 ) aud 
HCN(4-3) hues are practically the same. The frequeucies 
of the CO(3-2), HCN(4-3) aud CS(7-6) hues are close 
to each other aud beam parameters for them should be 
similar, too. However, the comparisou with the CO( 2 -l) 
is complicated by the siguificautly differeut frequeucies. 

The CO(3-2) hue is most probably saturated. Theu 
the HCN(4-3) hue should be saturated, too. The CS 
hue is weaker implyiug either a relatively lower optical 
depth (but close to uuity auyway) or a smaller size of the 
emissiou regiou. A simple modeliug usiug e.g. RADEX 
(vau der Tak et al. 2007) shows that the columu deu- 
sities of HCIN aud CS required to explaiu such optical 
depths are > 10^^ cm“^. The relative abuudauces of 
these molecules are < 10“^. Therefore the total gas col¬ 
umu deusity iu the clump is > 10^^ cm“^. Theu, us¬ 
iug the size estimated above we obtaiu the gas deusity 
n > 3 X 10^ cm“^. This estimate is cousisteut with the 
observatious of the HCN aud CS hues which require at 
least such deusities for excitatiou. 

The large HCN aud CS hue widths iudicate that the 
clump is gravitatioually uubouud. The viri al mass of 
this clump estimated iu the usual way (e.g. Ziucheuko 
et al.j|1994|) is Mvir ^ 30 Mq. At the same time, the 
uou-detectiou of dust emissiou implies au upper limit for 
mass orders of maguitudes lower. Therefore, the clump 
represeuts a trausieut eutity. 

It is worth meutiouiug the uou-detectiou of HCO^ aud 
N 2 H+ emissiou. Both molecul es cau be destroyed by 
dissociative recombiuatiou (e.g. Ziucheuko et al. 2009). 
Therefore their abseuce cau iudicate au euhauced louiza- 
tiou iu the deuse clump. 


5.2. Ionized gas 
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Figure 19. Lef t pa nel: the position-velocity map for the CO(3-2) emission as observed with the IRAM 30m telescope along the cut 
indicated in Fig. |l8| The contour levels are on logarithmic scale from 0.4 K to 31.5 K. The gap at Vlsr 24 kms“^ is caused by the 
emission at the reference offset position mentioned above. Right panel: the position-velocity map for the CO(3-2) emission as observed 
with the SMA along the same cut. The contour levels are on logarithmic scale from 0.5 to 11.5 Jybeam“^. 




Fig ure 20. Left panel: the position-velocity map for the HCO+(4-3) emission as observed with the SMA along the cut indicated in 
Fig. |l7| The contour levels are on logarithmic scal e fr om 0.4 to 7.6 Jybeam“^. Right panel: a part of the position-velocity map for the 
CO(3-2) emission measured by the SMA from Fig. |19| for the same intervals of the offset and velocity as in the left panel. 


In Fig. we plot the map of the 15 GHz continuum 
emission near the SMA1/SMA2 clumps (from the VLA 
archival data, the angular resolution is about 4 arcsec) 
which shows the distribution of the ionized gas in this 
area. This ionized component seems to be associated 
with the jet traced in the Fell emission. The GMRT 
map at 1280 MHz looks similar but there is a signifi¬ 
cant positional uncertainty in the GMRT data (Paper I). 
Properties of the continuum sou rce were estimated in Pa¬ 
per I and by |Ojha et al. (2011). The emission measure 


is EM ^ (1 — 2.5) X lO'^ pccm“^. Taking into account 
the observed size of the continuum source, the electron 
density is n ^ 3 x 10^ cm“^. 


6. SURROUNDINCS 

The surroundings of the SMA1/SMA2 clumps and high 
velocity outflow are traced in several molecular lines. 
One of the most informative is the N2H+(3-2) transition 
observed with both SMA and IRAM-30m telescope. The 


combined map of the N2H+(3-2) emission is presented 
in Fig. 18 It shows an absence of N 2 H+ in the hot core, 
in accoraance with our previous findings (Paper I). The 
overall morphology of the N 2 H+ emission suggests that 
it originates in an envelope around the central cores and 
the outflow lobes. 

The distribution of the SiO(5-4) emission is different 
(Fig. [^. It peaks near the SMAl core. There is also a 
feature m the area of the blue outflow lobe which may 
be associated with the N 2 H+ emission. The CS(7-6) 
distribution (right panel in Fig. 23) seems to be the most 
uniform one. It shows a rather smooth, almost spherical 
halo around the SMAl and SMA2 cores. There is no sign 
of the outflow influence on the CS distribution. One may 
suspect that the optical depth in the CS line is too high 
and we see only the outer regions of the core. However 
the SMAl and SMA2 clumps are well resolved in this 
map which makes such explanation less probable. Then, 
the total gas distribution is apparently not significantly 
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Figure 21. The CO(3-2), CO(2-l), HCN(4-3), CS(7-6), 

HCO+(4-3) andN2H+(3-2) spectra towards the position marked 
by cross in Fig. |18| 

affected, too. 


7. DISCUSSION 

The main goal of this study is the characterization of 
the outflow and probable accretion disk associated with 
the massive YSO in the S255IR clump. The data pre¬ 
sented in the previous sections shed new light on this 
system. 

Concerning the outflow, one of the main findings is 
that the SMA interferometric data alone give a rather 
distorted picture. They hint at two highly collimated 
parallel outflows with different centers of origin. How¬ 
ever, this impression is apparently caused by a signifi¬ 
cant flux loss in the interferometric measurements. The 
single-dish CO observations clearly show a single less col¬ 
limated outflow originating at the SMAl core. The CO 
emission retrieved by the SMA near the driving source 
originates apparently from the northern wall of the out¬ 



a (2000) 


Figure 22. The map of CO(3-2) emission in the velocity range 
-19...+1 kms“^ (contours) overlaid on image of CS(7-6) emission 
in the velocity range -13...+1 kms“^ (grey scale). The intensity 
units for the image are Jybeam“^ kms“^. The contour levels are 
(3, 5, 7, 9, 11, 13)X5 Jybeam“^ kms“^. The SMA beam is shown 
in the lower left corner. 


flow cavity. This means that this part of the wall is 
rather thin and bright. The absence of a noticeable CO 
emission at the opposite side of the wall implies a more 
diffuse distribution of the emission here, probably due to 
the density structure of the surrounding medium. 

The question arises of how common is this effect in 
the interferometric studies of outflows. We can easily 
imagine a situation when an interferometric image will 
show multiple outflows from a single driving source, while 
in fact there is a single wide-angle outflow. 

The obser ved P-V diagrams for the outflow 
(Figs. |19|20) most probably indicate periodic ejec¬ 
tions froni the driving source. Two events can be traced 
in the data. They apparently created jet knots at 
different distances from the central star. The older one 
is responsible for the extended CO outflow. The peaks 
of the CO emission coin cide with bright H 2 emission 
spots (Wang et al. 2011) which probably indicate bow 
shocks at the heads of the jets. The next ejection event 
created other bow shocks which are seen in particular in 
the Fell emission (Fig. [^. This later ejection entrains 
dense molecular gas observed in the wings of the HCO+, 
HCN and CS lines. It is also traced in the CO emission 
as ca n be s een in the CO posit ion-velocity diagram 
(Figs. |19|20). Another manifestation of this activ ity i s 
apparently the dense high velocity clump (Sect. |5 1[ ), 
which most probably represents dense gas at the head 
of bow shock. All the jet knots lay practically on a 
straight line. Therefore the orientation of the jet does 
not change with time significantly. The age of the first, 
ext ended outflow was estimated of about 7000 years 
by Wang et al. (2011). Assuming the same ejection 
velocities tor the two events and comparing distances of 
the bow shocks from the driving source we can conclude 
that the second ejection happened about 1000 years ago. 
Therefore the time interval between the ejection events 
is about 6000 years (the age estimates were done under 
the usual assumption of an inclination angle for the jet 
of 45°). 
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Figure 23. Left panel: the maps of the SiO(5-4) (left panel) and CS(7-6) (right panel) integrated line emission obtained by combining 
the SMA and 30m data (orange contours) overlaid on the continuum image at 0.8 m m. The maps of the CO(3-2) high velocity emission as 
observed with the IRAM 30m telescope (blue and red thick contours) from Fig. |18| are also shown. The beams for the IRAM-30m telescope 
and for the combined maps are shown in the lower left corner of each panel. 


Several entrainment mec hanisms are usually discussed 
for m olecular outflows (e.g. Arce et al.|2QQ7 Frank et ah 
'20141) . In our case the most probable on e is the jet bow- 

k ~ ~ - - - 


shock model (e.g. Raga & Cabrit 
gument in favor of this model is t 


1993|). The mam ar- 
le obvious presence of 


several bow shocks (traced in the H 2 and Fell emission) 
clearly associated with the high velocity molecular gas. 
This model is also consistent with the observed outflow 
morphology and kinematics. The shape of the P-V dia¬ 
gram differs from the frequently observed “Hubble law” 
and shows a range of velocities at the largest distances 
from the driving source. A similar shape was observed in 
some other outflows which are apparently driven by bow 


shocks and is expected in the theoretical models (e.g. Lee 
eTaLlpOPOl ). 


T'he jet is apparently launched from the accretion disk 
around the central massive young star. As mentioned 
above, there is a chain of water masers along the jet. For 
several water maser spots shown in Fig. p roper mo¬ 
tions have been measured (Goddi et al. 2007). The ve¬ 
locities (from ^ 10 kms“^ to ^ 30 kms“^ at the distance 
of 1.6 kpc) are mostly perpendicular to the jet axis with 
outward components in some cases. The maser veloc¬ 
ity pattern is consistent with clockwise rotation around 
the core center. Velocities of the masers closest to the 
center imply dynamical mass of the order of 20 Mq. 
The mass of the central star is estimated to be 27 Mq 
( jOjha et al.|| 2011 |) assuming the distance of 2.5 kpc. At 
the distance of 1.6 kpc the mass of the star would be 
23.5 ± 3.6 Mq. This estimate is based on the unpub¬ 
lished value of the star luminosity of ( 6.8 ± 2 . 8 ) x 10 ^ Lq 
( for the distance of 1.6 kpc). Now we reconsidered the 
star SEP by including data at millimeter wavelengths 
from Zinchenko et al.| ( |2009 ). This results in the lumi- 
nosity of (3.5 ± 0.3) x 1 Q^ ~L^ and mass of 20 ± 2 It 
is worth mentioning that Pashchenko et al. (2003) esti¬ 
mated the mass of the central star to be (6-y) Mq from 
their analysis of water masers variability. However, they 
employed a systemic velocity of 8.3 kms“^ which contra¬ 
dicts our findings which clearly show the systemic veloc¬ 
ity to be about 4.8 kms“^. In the following we will take 


the mass of 20 Mq as the most probable value. This 
is sufficient to explain velocities of the masers closest 
to the star. However maser spots located at larger dis¬ 
tances from the center have similar velocities and imply 
a much larger central mass of the order of 100 Mq. This 
value is inconsistent with all other estimates. The cor¬ 
responding mean density would be unrealistically high. 
This means that we probably need another explanation 
for the maser velocities. Rather speculative, it can be jet 
rotation. Rapid rotation o f protostellar jets is expected 
in theoretical models (e.g. Pnd ritz et al.|2QQ7t and was 
indeed observed in some cases (|Hacciotti et al.| 2002). 

The water masers are thought to be excited by inter- 
action of the jet with the surrounding medium. The fact 
that they are moving in one direction on each side of 
the central star implies some asymmetry in the system. 
There may be some misalignment between the jet and 
the rotation axis of the material in the outer parts of 
the clump (apparently, in the axially-symmetric case the 
observed maser movements would be symmetric relative 
the axis). It can be that in the inner part the orientation 
of the disk plane is different. However, the difference 
cannot be large. We see also other indications of the jet 
interaction with the surrounding medium, in particular, 
increased line widths along the jet. 

Fig. shows that the water maser condensations are 
shifted from the jet axis in direction of their proper mo¬ 
tion. This projected shift is of the order of 200 AU. At 
the velocities mentioned above, the maser condensations 
would travel this distance in about 30-100 years. From 
this consideration it looks probable that the event cre¬ 
ated these masers happened quite recently. 

The disk is apparently strongly fragmented as follows 
from the derived beam filling factor (^ 0.2) for CH 3 CN 
and CH 3 OH. The mean gas density in the disk is about 
6 X 10^ cm“^. The density of the fragments should be 
much higher. Observations of vibrationally excited HCN 
indicate densities > 10 ^^ cm“^ which is consistent with 
this picture. 

The derived mass of the hot gas (^0.3 Mq) is by an 
order of magnitude insufficient to gravitationally bound 
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the core. The virial mass estimated in the usual way 
(e.g. |Zinchenko et al.||1994|) from the line width is about 
7 Mq. I'here can be also some amount of cold mate¬ 
rial in the disk, as follows from the DCN observations. 
However, the total mass of the disk is apparently much 
lower than mass of the central star which is estimated 
of about 20 Mq (see above). Since the disk mass is sig¬ 
nificantly lower than the star mass, the disk is probably 
not self-gravitating and we may expect Keplerian rota¬ 
tion. Unfortunately the achieved angular resolution is 
not sufficient to measure reliably the disk rotation curve. 

For a central mass of 20 Mq the inclination angle of the 
disk (derived from comparison of the observed rotation 
line-of-sight velocity with the expected velocity for this 
mass) should be small, about 25° (Fig. [^, i.e. the disk 
is seen almost face-on. This is consistent with our maps 
which do not show a significant elongation of this object. 
On the other hand this conclusion looks somewhat sus¬ 
picious because of the well-collimated appearance of the 
molecular outflow. Observations at a higher angular res¬ 
olution are needed to clarify the structure and kinematics 
of this object. 

The temperature of the hot gas derived from the 
CH 3 CN, CH 3 OH and SO 2 observations is about 130- 
180 K. A similar estimate for the temperature of the 
hot component was obtained by us many years ago from 
the analysis of the IRAS data (Zinchenko et al. 1990). 
The HNCO rotational temperature is somewhat higher, 
about 320 K. However, as mentioned above, this value 
probably represents an upper limit to the excitation tem¬ 
perature. In addition, the HNCO excitation can be in¬ 
fluenced by the FIR radiation. The HNCO excitation in 
massive cores was discussed by Zinchenko et al. (2000). 
That analysis shows that in the present case colhsional 
excitation can be effective. The critical density for exci¬ 
tation of the K_i = 3 ladder is of the order of 10^^ cm“^. 
Such density is quite possible in the SMAl core as fol¬ 
lows from the estimates above. However the radiative 
excitation by FIR emission via the 6 -type transitions 
cannot be excluded, too. With the column density of 
A^(H 2 ) ^ 3 X 10^^ cm“^ the dust will be optically thick 
at these frequencies creating a sufficiently strong radia¬ 
tion field. The HNCO excitation especially in the higher 
K-i ladders may reflect the effective temperature of this 
field. 

The highest temperature of the observed molecular 
material is about 300 K. Taking the source luminosity 
to be about 3.5 x 10^ Lq (see above) this value corre¬ 
sponds to the dust equilibrium temperature at the dis¬ 
tance of 200-300 AU from the star. Since we do not see 
warmer molecular gas, no dust should be present within 
this radius in case of no shielding. However this radius 
is comparable to the observed size of the core and we 
do not see any central hole in the molecular distribution. 
It is then possible that shielded molecular clumps are 
distributed at smaller radii. 

In this respect an interesting feature of our observa¬ 
tions is the apparently relatively cold and rather massive 
clump with a strong DCN emission. It may be gravita¬ 
tionally bound and the mass of the clump is sufficiently 
high to consider it as a possible low-mass protostar. 

Concerning the surroundings of the SMAI clump, we 
see that the outflow strongly affects the chemical compo¬ 
sition of the medium. The N 2 H+ molecules are destroyed 


along the outflow. The SiO distribution seems to be also 
significantly affected. At the same time there is no no¬ 
ticeable influence on the CS distribution. The density 
structure is probab ly not affected, too. As mentioned 
above (Sect. 4.2|4.4) the nearby clumps are probably in¬ 
fluenced by shocks which are apparently associated with 
this outflow. 


8 . CONCLUSIONS 

We presented the results of our observations of the 
S255IR area with the SMA at 1.3 mm in the very ex¬ 
tended configuration and at 0.8 mm in the compact con¬ 
figuration as well as with the IRAM-30m at 0.8 mm. 
The best achieved angular resolution is about 0.4 arcsec. 
The dust continuum emission and several tens of molecu¬ 
lar spectral lines are observed. The majority of the lines 
is detected only towards the S255IR-SMAI clump. In 
summary, our main findings are the following: 

1. The S255IR-SMAI clump represents apparently a 
rotating structure (probably a disk) around the young 
massive star. The achieved angular resolution is still in¬ 
sufficient to establish the character (Keplerian or non- 
Keplerian) of the rotation. The temperature of the 
molecular gas reaches I30-I80 K. The size of the clump 
is about 500 AU. It is apparently strongly fragmented 
as follows from the derived small (^ 0 . 2 ) beam filling 
factors for various molecules. The mean gas density is 
about 6 X 10^ cm“^. The density of the fragments should 
be much higher which is confirmed by observations of 
HNCO and vibrationally excited HCN. The mass of the 
hot gas is ^ 0.3 Mq and the total mass of the clump 
is significantly lower than the mass of the central star 
(about 20 Mq). The inclination angle of the disk (de¬ 
rived from comparison of the observed rotation line-of- 
sight velocity with the expected rotation velocity for this 
mass) should be small, about 25°. 

2 . We detected a strong DCN(3-2) emission near the 
center of the SMAI clump. Most probably it indicates 
the presence of a rather large amount (> I Mq) of cold 
(< 80 K) material. This cold clump can be gravitation¬ 
ally bound. 

3. High velocity emission is observed in the CO line 
as well as in lines of high density tracers HCN, HCO+, 
CS and other molecules. The CO outflow is much more 
extended than that observed in the lines of high density 
tracers. Its morphology obtained from combination of 
the SMA and IRAM-30m data is significantly different 
from that derived from the SMA data alone. The CO 
emission detected with the SMA traces only one bound¬ 
ary of the outflow and leads to a rather distorted picture 
of the outflow structure. The velocity of the CO outflow 
reaches ^ 60 kms“^. 

4. The outflow is most probably driven by the jet bow 
shock mechanism. The available data indicate at least 
two major ejection events with a time interval of sev¬ 
eral thousand years between them. The direction of the 
ejections does not change with time. The high velocity 
emission in the lines of high density tracers is associated 
with the peaks of the Fell emission related to the bow 
shocks caused by rather recent ejections from the SMAI. 

We detected a dense high velocity clump associated 
apparently with one of the bow shocks. It shows a strong 
emission in the HCN(4-3) and CS(7-6) lines. At the 
same time there is no detectable HCO+(4-3) or N 2 H+( 3 - 
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2) emission which can be probably explained by enhanced 
ionization. 

5. The proper motions of the water masers excited 
along the jet imply some misalignment of the jet with 
the rotation axis of the material in the outer parts of the 
clump. It can be that the orientation of the disk in the 
inner and outer parts is somewhat different. However, 
the difference cannot be large. 

6. The outflow strongly affects the chemical composi¬ 
tion of the surrounding medium. The N 2 H+ molecules 
are destroyed along the outflow. The SiO distribution 
seems to be also significantly affected. At the same time 
there is no sign of the outflow influence on the CS dis¬ 
tribution. The total gas distribution is apparently not 
significantly affected, too. 
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